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From general arguments, that are valid for spin models with sufficiently short-range interactions, 
we derive strong constraints on the excitation spectrum across a continuous phase transition at zero 
temperature between a magnetic and a dimerized phase, that breaks the translational symmetry. 
From the different symmetries of the two phases, it is possible to predict, at the quantum critical 
point, a branch of gapless excitations, not described by standard semi-classical approaches. By 
using these arguments, supported by intensive numerical calculations, we obtain a rather convincing 
evidence in favor of a first-order transition from the ferromagnetic to the dimerized phase in the 
two-dimensional spin-half model with four-spin ring-exchange interaction, recently introduced by 
A.W. Sandvik et al. [Phys. Rev. Lett. 89, 247201 (2002)]. 
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I. INTRODUCTION 

In the last years, a large amount of work has been 
done to clarify the properties of unconventional magnetic 
systems. Indeed, the presence of frustrated interactions, 
generated either by the geometry of the lattice or by the 
competing interactions, can give rise to many anoma- 
lous low-energy properties^ One of the most exciting is 
the possibility to obtain phases with no magnetic order 
even at zero temperature^ The actual interest in this 
subject has been recently renewed by the discovery of 
several materials^^ that do not show any sign of mag- 
netic order down to very low temperatures. From a the- 
oretical point of view, a considerable progress has been 
done in elucidating the possible unconventional states 
and, apart from magnetic phases, it is now well accepted 
that many systems exhibit a large variety of disordered 
ground states. One class of such paramagnetic states is 
given by the so-called valence-bond solids, where pairs of 
nearest- neighbor spins form a singlet, leading to an or- 
dered pattern of valence bonds. These states explicitly 
break the translational symmetry, implying a degenerate 
ground state, and represent the two-dimensional exten- 
sion of the celebrated dimerized state, rather well estab- 
lished in quasi-one-dimensional materials. 6 Another class 
of disordered phases is given by the so-called spin liquids, 
that do not break any local symmetry. It is now rather 
widely accepted that also spin liquids possess some de- 
generacy, or even a gapless spectrum, and have some kind 
of order, called topological order, related to non-local op- 
eratorsi^ 

Much recent research has focused on the way to de- 
scribe transitions between different quantum phases. In- 
deed, from the Landau's arguments valid for classical 
critical phenomena, if one exclude a delicate fine tuning 
of the parameters, there is no reason to have a continuous 
transition between two phases with different symmetries 
of the order parameters, like for instance between a mag- 



netic system and a valence-bond solid, and a first-order 
transition would be the most natural scenario. However, 
very recently^S*^ these arguments have been questioned 
and it has been argued that the situation can be more 
complex and richer in quantum systems, where the Lan- 
dau paradigm could be violated, giving rise to an appeal- 
ing way to escape a first-order transition. Therefore, it 
has been suggested that a quantum phase transition is 
intrinsically different from a classical one. In particular, 
when considering quantum phases, the natural descrip- 
tion is no longer by using order parameters but instead 
by considering new degrees of freedom carrying fractional 
quantum numbers, that emerge at the critical point. 

The arguments given in Refs. ©and US based on field- 
theory approaches, are rather convincing but it is natural 
to ask if microscopic spin models, defined on a lattice, 
could present the same low-energy critical behavior. In 
this sense, however, we face a technical difficulty: a good 
diagnosis would require to study large enough clusters, 
but this is not possible since quantum Monte Carlo sim- 
ulations of frustrated antiferromagnets are plagued with 
a very severe minus sign problem. In this respect, it is im- 
portant to find out microscopic models that do not suffer 
from this problem but, at the same time, present phase 
transitions between quantum phases. Recently, Sandvik 
and collaborators^ introduced a spin-half model on the 
two-dimensional square lattice in the presence of a par- 
ticular ring-exchange interaction: 

ti = TixY + HRjngj (1) 

where 

n mng = -k {s+s^s+sr + srs+StS+). 

[i,j,k,l] 

Here denotes a pair of nearest neighbor sites and 

[i,j,k,l] indicates the sites on the corners of a plaque- 
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tte. In particular, Wrung is a ring-exchange interactions 
that acts as a cyclic permutation on the spins of a given 
plaquette whenever they are in a S z = configuration 
with parallel spins on opposite corners^ The case with 
J > and K > is considered, in order to avoid the sign 
problem. Notice that H is invariant under global U(l) 
transformations — ► e ie . 

The Hamiltonian Q has been investigated by using 
quantum Monte Carloii*^ and it has been shown that, 
at zero temperature and zero magnetic field, there is 
a transition between a ferromagnetic phase, stable for 
small values of Kj J, and a dimerized one, characterized 
by columnar order at Q — (tt, 0) and Q — (0, n), sta- 
ble for intermediate ring-exchange couplings. Finally, for 
large K/J a checkerboard phase is also expectedj 1 ^* 1 ^ 
Note that, in our notation, J is twice time bigger than 
the one used in Refs. [ll] andEl Once scaled in our 
units, their best estimate for the transition point from 
the magnetic to the dimerized phase is K c /J ~ 3.957. 14 
In the original paper^l it has been shown that, by in- 
creasing the strength of the ring exchange, both the spin 
stiffness p s and the magnetization M decrease and even- 
tually vanish, with a corresponding insurgence of a finite 
dimer order parameter. Hence, the Hamiltonian Q has 
been suggested as a prototype model to test the ideas 
of Refs. El and I lfj However, subsequent and more accu- 
rate calculations have raised some doubt on the previous 
ones, suggesting instead a weak first-order transition. 15 
The actual nature of the transition is therefore still con- 
troversial. 

In this paper, we propose a different point of view 
to investigate the properties of a phase transition be- 
tween a magnetic and a dimerized phase of a generic 
spin Hamiltonian with sufficiently short-range interac- 
tions. Indeed, instead of calculating the order parameters 
and the stiffness, we derive strong constraints on the ex- 
citation spectrum across a continuous phase transition 
between a magnetically ordered and a dimerized phase, 
breaking translational symmetry. These constraints can 
be verified by using numerical techniques, such as the 
Green's function Monte Carlo method^ In our opinion 
this way to proceed represents the first attempt to clar- 
ify the evolution of the excitation spectrum and can give 
important insight into the nature of the transition point 
also in other contexts. 

The paper is organized as follow: in Sec. [U] we develop 
the general theoretical formalism, in Sec. 11111 we present 
our numerical results, and in Sec. II VI we draw our con- 
clusions. 



II. THEORETICAL FORMALISM 

In this section, we describe in detail the implications 
of a continuous phase transition from a magnetically or- 
dered phase to a dimerized one on the excitation spec- 
trum near the quantum critical point. In the following, 
we will consider an Hamiltonian with a generic transla- 



tional invariant two-spin coupling and SU(2) symmetry, 
and at the end of the section we will also discuss the cases 
with U(l) symmetry or multi-spin interactions. 

Let us consider the Heisenberg Hamiltonian on a two- 
dimensional lattice with N = L x L sites: 

H = ^2Ji,jSi-Sj, (2) 

where Ji j is the coupling between the spins at site Ri 
and Rj and Si = (Sf,Sf, Sf) is the spin operator at the 
site Ri. For a translational invariant magnetic coupling 
that depends only upon \Ri —Rj\, i.e., Jij — J\i-j\ = J r , 
we can perform the Fourier transform and easily obtain: 

1~L = ^ JqSq • S—qj (3) 
9 

where S q = 1/VN £ . e lqR i Sj and J q = £ r J r de- 
fine the Fourier transform of the local spin operator and 
the magnetic coupling, respectively. In the following, we 
will consider interactions that are sufficiently short range, 
that is: 

j=^J2\ J i,i\\Ri-Ri\ 2 <™- (4) 

For the SU(2) case, the eigenstates can be classified ac- 
cording to the total spin S and the ground state, here 
denoted by \^q ), is generally a singlet with momentum 
Qo. By means of the so-called Feynman construction, 
originally introduced to describe the excitations of the 
liquid Helium, 17 we can easily define a variational state 
for the triplet excited state with momentum q with re- 
spect to the ground state: 

|* Qn+g ) = S%\* Qo ), (5) 

S™ being one of the three components of the total spin S q 
(a = x, y, z). Then, it is easy to calculate the variational 
energy E Qo+q of |*Q 0+g > 

A(Q + q)= E Qa+q - E Qa = ^ q ga) , (6) 

where (...) stands for the average value over |^q ) and 
the function is given by the double commutator 

K = 1([ s -i>[h>^}})- ( 7 ) 

By using standard commutation relations for the spin 
operators and by summing over all the spin components 
a, we obtain 

E F , = 2^ £ [ Jfc +« + Jk -v ~ 2 Jfc l ^ ■ S-x) ■ ( 8 ) 

a k 
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After performing the inverse Fourier transform, we have 

E K = Jf E JiA^[q{Ri - R,)} - l}(Si ■ S 3 ) < 

a i,j 

£fY^\Jv\\Ri-Rj\ 2 \0i-Si)\< JS(S +!)<?, (9) 

i,j 

where we used that \{Si ■ Sj) \ < S(S + 1) and the short- 
range condition . Notice that the previous bound does 
not rely upon a small-q expansion, but it is instead valid 
for all momenta. Therefore, we arrive to the important 
result that: 

where S(q) = (S" q Sq) is the static magnetic structure 
factor, that, in the SU(2) case, does not depend upon a. 
Since \^Q 0+q ) is a variational state, its average energy 
has to be higher than the lowest exact triplet excitation 
with given momentum Qq + q. We conclude, therefore, 
that the inequality (|1L)|> will also hold for the exact triplet 
excitation energy with momentum Qq + q. 

In the following we will argue that the rigorous bound 
on the excitation spectrum IjlUI) quite naturally implies 
the existence of a new branch of gapless excitations close 
to an hypothetical second-order phase transition from a 
magnetic phase to a dimerized one. We will focus on the 
square lattice case and an ordered phase with staggered 
magnetization at Q = (tt,tt), but similar results can be 
also found in the case of different lattice structures or 
spiral magnetic orders. For the sake of simplicity, we can 
imagine that J q depends upon one parameter g, propor- 
tional to the frustration of model: for instance, in the 
J\ — J 2 model J q = 2 (cos (/j; + cos^j,) -I- 4g cos q x cos q y , 
with g = J21 1 J\- Suppose that, by increasing the frustrat- 
ing parameter <?, the system has a continuous transition 
from a magnetically ordered phase to a dimerized one. In 
this latter, the lowest triplet excitation remains gapped 
in the thermodynamic limit and the ground state is four- 
fold degenerate, with four different momenta q — (0,0), 
(0, 7r), (tt, 0), and (tt,tt). Again, the following arguments 
will also hold for any other dimer patterns, having ground 
states with momenta q 7^ (0,0) and q ^ (tt, tt). By ap- 
plying the Feynman construction with q = (tt, tt) to the 
state with q = (0, tt) and by using the inequality ljTU|) . we 
obtain that 

o(TT, TT) 

In the dimerized region, i.e., for g > g c , the spectrum 
is gapped, implying that S(tt, tt) is finite, and the con- 
straint l|llf> does not give us any useful information on the 
excitations. However, if the transition is continuous, the 
four-fold degeneracy will hold down to the critical point 
g = g c , where S(ir,n) diverges, because of the incipient 
magnetic phase, stable for g < g c . Indeed, at the critical 
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FIG. 1: (Color online) Comparison on a 6 x 6 lattice between 
Lanczos and GFMC for the lowest energy gap. Upper panel: 
A + (-7r, 0) (see text) [in unit of (K + J)], as a function of K/J. 
Lower panel: Asp(n,0) (see text) [in unit of (K + J)], as 
a function of K/J. Inset: direct comparison of the dynami- 
cal correlation function for the dimer operator calculated by 
GFMC approach and by Lanczos. 

point the spin-spin correlations are expected to decay as 
-(!+»?), leading to S{tt, tt) - L 1 ^ . The actual value 
of 77 can be found by considering the corresponding clas- 
sical model in three dimensions. In the Heisenberg model 
with nearest-neighbor interactions, we have that 77 ~ 0, 
implying a rather strong divergence of the static struc- 
ture factor at the critical point. However, this divergence 
can be much weaker, i.e., rj ~ 0.5, in the case where no 
free topological singularities are allowed but only pairs of 
them are present li&iS We point out here that our conclu- 
sions remain valid as long as rj < 1, the structure factor 
being logarithmically divergent in the case of i] = 1. 

Therefore, because the operator S% used in the Feyn- 
man construction carries S = 1, we immediately arrive 
to the conclusion that Eq. (|ll|l implies the existence of 
a branch of excitations, with S = 1 and q = (jr, 0), that 
becomes gapless at g = g c . Analogously, by considering 
the ground state with momentum q = (tt, 0), we obtain 
a branch of excitations with S = 1 and q = (0,7r). This 
outcome contrasts the common understanding of the low- 
energy properties of magnetic systems with 0(3) symme- 
try, which contains only two branches of gapless excita- 
tions at q = (0,0) and q — (tt,tt), while the triplet exci- 
tations at q = (0, tt) and q — (0, tt) are always gapped^ 
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FIG. 2: (Color online) The dimer gap Asp(7r, 0) [in unit of 
[K + J)] obtained by GFMC as a function of the size N of 
the cluster for different values of K/J. The lines are linear 
fits of the points. Inset: the extrapolated value in the ther- 
modynamic limit as a function of K/J. 

Therefore, one is left with the following three possibili- 
ties: the first one is to have an 0(3) critical point with 
more than two gapless modes, not described by standard 
theories of magnetic phase transitions. The second one is 
that the dimerized phase is separated from the magnetic 
phase by a first-order transition: in that case the spec- 
trum can evolve discontinuously and the triplet states at 
q = (0, 7t) and q = (0, tt) can remain gapped. The third 
possibility is that there is still a continuous transition but 
the non-magnetic phase does not break the translational 
invariance and, therefore, has no dimer order. Then, the 
dimerized phase could be eventually stabilized through a 
further (continuous) transition. 

We stress that a continuous transition between a mag- 
netically ordered phase and a dimerized one is instead 
possible in models that explicitly break the translational 
invariance, like for instance the coupled ladder system^Si 
where due to the folding of the Brillouin zone q — (tt, 0) 
[or q = (0, 7r)] is equivalent to q — (0, 0). 

Now we would like to discuss what happens for systems 
with different kinds of interactions. In the case of multi- 
spin interactions, it is easy to prove that we can arrive to 
similar expressions and the only difference is the prefac- 
tor, due to a different value of the double commutator J7J). 
A little more care must be paid to a model with U(l) sym- 
metry. Here, the total spin S does not commute with the 
Hamiltonian, and only S z is a good quantum number. 
Furthermore, the states with S z = can be also classi- 
fied according the discrete particle-hole symmetry (PH) 
S z — > — S z . In the weakly frustrated limit, without loss of 
generality, we can consider the ferromagnetic XY model, 
that has only one gapless branch at q = (0, 0). On a finite 
cluster, the ground state has S z — and can be taken 
to be even under the PH symmetry, i.e., with PH = +1. 
In this case, we can construct two different variational 
states for the low-energy excitations. The first one, with 
S z = 1, is given by |* Qo+9 ) = S+\* Qo ), S+ being the 



Fourier transform of = SJ + iSj, and has 

Instead, the second possibility, with S z — and PH = 
-1, is given by \^ Qo+q ) = S^ Qo ) and has 

A*(Qo + 2) = ^. (13) 

In the dimerized phase the ground state is four-fold de- 
generate with four different momenta q = (0,0), (0,7r), 
(tt, 0), and (tt, tt), all of them having S z — and PH = 
+1 [this is equivalent to have all singlets in the SU(2) 
case]. Then, in close relation with the SU(2) case, the 
Feynman construction l|12() with q — ► (0, 0) applied to 
the state with momentum q — (tt, 0) implies the exis- 
tence of a gapless branch of excitations at q = (tt, 0) 
with S z 7^ 0. Indeed, in this case F+ is finite for \q\ — > 
and the structure factor in the denominator diverges at 
the critical point, as (SZ q S q ) oc l/q 1 ^ 71 . Hence, we ar- 
rive to the same conclusions as before: excluding the 
possibility of having more than one gapless branch in 
a model with 0(2) symmetry, in contrast with semi- 
classical approaches^ the transition must be either first 
order or second order to a non-dimerized phase. Fi- 
nally, it should be noted that, under an additional as- 
sumption, also Eq. (|13f) gives a strict bound to the spec- 
trum. Indeed, F z oc q 2 for |g| — > and, whenever at 
the critical point (SZ q S z ) does not vanish too fast, i.e., 
(Sl q S z ) oc q 2 - a with a > 0, also A z (q) for \q\ -> 0. 

III. NUMERICAL RESULTS 

In this section, we present our numerical results for the 
Hamiltonian Q), by using the Green's function Monte 
Carlo (GFMC) technique with a finite and fixed popu- 
lation of walkers4& In order to minimize the statistical 
fluctuations, we implemented the importance sampling, 
defined by the variational wave function^ 

\*v)=V t J\F), (14) 

where \F) is the ferromagnetic state with the magnetiza- 
tion along the x direction, i.e., \F) = Ui(\ |)j + | J.),), V z 
is the projector onto the subspace with S z = 0, and J is 
a spin Jastrow factor 

^ = ex p j^E^A 25 jj> (15) 

with the two-body potential Vij that depends only upon 
the distance \Ri — Rj\ and can be optimized by using the 
minimization technique described in Ref. l23l 

Since the Hamiltonian does not suffer from the sign 
problem, it is possible to sample the exact ground-state 
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FIG. 3: (Color online) The spin-wave spectrum [in unit of 
(A' + ,/)] obtained by GFMC by using the operator O q = S%. 
The linear sizes of the cluster are L — 16 and 18, and the 
lines are guides to the eye. The vector q is taken along the 
(1,0) direction. Upper inset: the behavior of the finite-size 
gap at the minimum momentum, i.e., q m in = (2tt/L, 0) as a 
function of L for different values of K/J. Lower inset: the 
thermodynamic limit of the gap at q — (0, 0) as a function of 
K/J. 



wave function <&o(x) — (x|$o) (being \x) a generic spin 
configuration), or, with importance sampling, the posi- 
tive quantity <J> (a;)<I>y {%) = (x\®o)(x\®v)- Moreover, it 
is possible to evaluate the imaginary time evolution of 
the following dynamical correlation function: 



G{q,r) 



(§v\O- q e- rH O q \* ) 



(16) 



where O q is a given operator with a definite momentum 
q. This technique has been already successfully used by 
one of us to compute the excitation spectrum of a quasi- 
one-dimensional magnetic system, 24 It is easy to show 
that, for large imaginary time r, G(q,r) behaves like: 



G(q, t) oc e 



(17) 



where A q = E q —Eq is the energy gap between the ground 
state and the first excited state |$ ) such that 



($ g |O 9 |$ > ^0. 



(18) 



Therefore, from the fit of the large-r behavior of the dy- 
namical correlation function, it is possible to extract the 
lowest gap with a given momentum q and the symmetry 
properties imposed by the condition (TSJ. 

Whenever the operator O q is diagonal in the spin con- 
figuration \x), like, for instance the case of S z , the ex- 
pression (|16fl can be easily calculated by using the stan- 
dard forward walking technique: 16 On the other hand, for 



non-diagonal operators, i.e., for S+, a more involved cal- 
culation is needed, since the operator breaks the Markov 
chain of the Monte Carlo simulation. In order to avoid 
this extra numerical effort, we have calculated more effi- 
ciently A + (q) by sampling the ground state with S z = 1 
and by calculating its excitation spectrum, for instance 
by using S q . Then, A + (q) can be found from: 



A+(q) = E q 
= (E, 



i 



-E, 



S z =l\ 



■ E; 



s*=i 



E, 



s z =o 



,(19) 



where the term in the first bracket can be calculated by 
fitting G(q, r) and the last two terms are just the ground- 
state energies of the two sectors with S z — and S z = 1 , 
easily calculated by GFMC. 

In order to verify the accuracy of GFMC for the eval- 
uation of the gap, we can perform calculations on a 6 x 6 
lattice, where the exact results are available by the Lanc- 
zos method. In particular, we consider the dimer opera- 
tor Oq - O s q p = e iqRi S?S? +x with q = (tt, 0). 
In this case, the operator O q does not change the PH 
quantum number and we have access to the lowest gap 
with the same PH of the ground state. In Fig. ^ we 
compare, for different values of the ring-exchange cou- 
pling K/J, the exact energy gap at q = (tt, 0) and even 
under PH with the one extracted from the large-r be- 
havior of G(q,r), and denoted by Asp(ir,0). Moreover, 
as a further check of the statistical approach, we can also 
compute exactly the dynamical correlation function and 
make the comparison with the GFMC calculations (see 
the inset of Fig. 0. It should be stressed that, for our 
purpose, we are only interested in the lowest energy gap, 
and, therefore, we fit only the large-r part of G(q,r). 
This can be achieved by using only two or three expo- 
nents at mostiS^ In the same figure, we also report the 
comparison for A + (7T,0), extracted from the decomposi- 
tion of Eq. (|19(l : also in this case the agreement is excel- 
lent. 

The lowest dimer gap Asp(tt, 0) provides us the in- 
formations on the appearance of a dimerized phase with 
spin-Pcicrls order. In Fig. we report the size-scaling 
calculations of Asp(7r, 0) as a function of the inverse of 
the size N, for different values of the ring exchange. The 
dimer gap is clearly finite in the magnetically ordered 
phase, with a decreasing behavior by approaching the 
transition. The thermodynamic limit indicates a vanish- 
ing value of Asp(7t, 0) for K/J > 3.95, in close agreement 
with Ref.EH 

Let us now turn to the spin- wave operators S q and S z , 
in order to apply the arguments developed in the previous 
section and assess the nature of the transition. First of 
all, it should be noticed that, in the XY model, S* acts 
on the ground state by creating a standard spin-wave 
excitation. This can be easily seen by considering the 
Holstein-Primakoff representation at the leading order in 
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FIG. 4: (Color online) The gap A+(tt,0) [in unit of (K + 
J)] calculated by GFMC as a function of the size N of the 
cluster for different values of K/ J. The lines are linear fits of 
the points. Inset: the thermodynamic value of the gap as a 
function of K/J. 

the 1/S expansion: 

S% = S-afa 

s i = — *yj ( a i - 

After standard calculations, we arrive to the result that 
S»|0) ~ aJ|o) 
S*\0) ~ a t| ), 

where |0) is the spin- wave ground state and a q is the cre- 
ation operator of the elementary excitation, defined by 
the Bogoliubov transformation. Therefore, the states of 
the gapless branch have non-zero overlap with the (nor- 
malized) states generated either by S% or by Sg, allowing 
us to assess directly the spin-wave spectrum. As shown 
in Fig (J3J), for small values of the ring-exchange coupling, 
the excitation spectrum remains gapless for \q\ — > and 
the opening of a gap, for K/J ~ 3.95, signals the tran- 
sition to the disordered phase. Moreover, in the ordered 
phase, the spectrum is linear for small momenta, i.e., 
A 2 (g) ~ c\q\ (see Fig. |3J), the constant c defining the 
spin-wave velocity. From our numerical results, it comes 
out that, as assumed in semi-classical approaches, the 
spin velocity c is not strongly renormalized by frustra- 
tion and remains finite up to the transition. 

By using the same spin-wave operators S£ and at 
q = (jr, 0) and q = (0, 7r), we can assess the nature of the 
transition between the ferromagnet and the dimer state. 
In particular, according to Eq. (|12l) . a continuous phase 
transition implies a vanishing A + (g) for q — (tt, 0) and 
q — (0, 7r). In Fig. (gj, we report the results for the gap 
A + (7r, 0) for different ratios K/J and various sizes of the 
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FIG. 5: (Color online) The gap A z (tt,0) [in unit of (K + 
J)] calculated by GFMC as a function of the size N of the 
cluster for different values of K/J. The lines are linear fits of 
the points. Inset: the thermodynamic value of the gap as a 
function of K/J. 

cluster. From our results, it turns out rather clearly that 
A + (7r, 0) remains finite in the thermodynamic limit for 
all the values of the ring-exchange couplings, also very 
close to the transition point. Of course, we cannot rule 
out the existence of a very tiny region, near to the tran- 
sition point, in which this gap would drop down to zero. 
Another possibility for missing the vanishing behavior of 
the gap could be due to an anomalous size scaling, that 
changes its behavior beyond a characteristic length scale 
L c much bigger than the values of L considered here. 
However, from our GFMC calculations, we do not have 
any sizable sign that could indicate the existence of gap- 
less excitation at q — (tt, 0) or q = (0, tt), and even at the 
transition point K c / J ~ 3.957 the extrapolated value of 
A + (7r, 0) remains clearly finite. Since a clear dimer order 
appears as soon as we enter the disordered regionjii we 
can safely rule out the possibility to have a spin-liquid 
phase in between the ferromagnetic and the dimerized 
phases. Therefore, our results on A + (7r,0) strongly sug- 
gest that a first-order transition is the more likely sce- 
nario in this 0(2) model with ring-exchange interactions. 
Similarly, also the gap A z (g) at q = (n,0) and q = (0,7r) 
does not show any evidence of softening and remains fi- 
nite at the transition, see Fig. [5] 

IV. CONCLUSIONS 

In conclusion, we investigated the excitation spectrum 
of a spin model close to the transition from a magnetic 
phase to a dimerized one, which is particularly impor- 
tant in view of the theoretical predictions of Senthil and 
coworkers. In particular, we considered an XY model 
in presence of a frustrating ring-exchange interaction, a 
model that certainly leads to a magnetic-dimer transi- 
tion, as recently shown by Sandvik and collaboratorsiiiii^ 

The central part of our work is the derivation of a gen- 
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eral constraint for the excitation spectrum for a broad 
class of short-range spin Hamiltonians: whenever the 
transition is second-order a new branch of gapless ex- 
citations must exists, besides the ones expected from 
standard semi-classical approaches. In particular, on the 
square lattice with a typical four-fold dimcr ground state, 
these modes are spin-wave-like excitations at momenta 
q = (0, 7r) and q — (tt, 0). On the contrary, by using 
a numerically exact technique, we did not find any ev- 
idence of these gapless modes. Therefore, our results 
rule out both conventional and unconventional continu- 
ous transition, implying that a first-order transition rep- 
resents the most likely scenario, in agreement with recent 



calculations. 15 Although we cannot exclude anomalous 
finite-size effects, that should be particularly important 
for correlation functions (see Ref. fill) , we expect that the 
excitation spectrum does not suffer from these problems, 
depending only on total energy differences. 
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